Exact Z 2 scaling of pair production in the high-energy limit of 

heavy-ion collisions 
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The two-center Dirac equation for an electron in the ex- 
ternal electromagnetic field of two colliding heavy ions in the 
limit in which the ions are moving at the speed of light is ex- 
actly solved and nonperturbative amplitudes for free electron- 
positron pair production are obtained. We find the condition 
for the applicability of this solution for large but finite colli- 
sion energy, and use it to explain recent experimental results. 
The observed scaling of positron yields as the square of the 
projectile and target charges is a result of an exact cancel- 
lation of a nonperturbative charge dependence and holds as 
well for large coupling. Other observables would be sensitive 
to nonperturbative phases. 

PACS number: 34.50.-s, 25.75.-q, 11.80.-m, 12.20.-m 

There is a very small number of problems in physics 
that yield to an exact solution. Remarkably, electro- 
magnetic production of free electron-positron pairs in 
the high-energy limit of peripheral, heavy-ion collisions 
can be described by a two-center, time-dependent Dirac 
equation which can be solved exactly and in closed form 
0. In this Letter, we study this exact solution and dis- 
cuss its implications for recent experiments performed at 
CERN's SPS 0, and possible future experiments at new 
facilities such as BNL's RHIC and CERN's LHC. (See 
Refs. [||JJ] for relevant reviews of this field.) 

Perturbative calculations have been held suspect at 
high energies and for the heaviest projectiles, e.g. because 
the coupling constant is not small (Za ~ 0.6) It 
is therefore surprising that positron yields observed from 
pair production in peripheral collisions of Pb 82+ ions at 
33-TeV on a Au target scale as Z\Zp, and that the ob- 
served positron-momentum distributions display an over- 
all good agreement with the leading-order perturbation- 
theory calculations Q . The exact nonperturbative solu- 
tion presented here explains these effects and is also con- 
sistent with the observed enhancement of the positron 
yields at small positron momentum. For future exper- 
iments, we indicate what observables would show com- 
plete agreement with second-order perturbation theory, 
and what other observables should be measured in order 
to detect nonperturbative effects. 

The relativistic scattering problem of an electron in the 
external field of two point-like charges (ions) , moving on 
parallel, straight-line trajectories in opposite directions 
at speeds which approach the speed of light, and at an 
impact parameter 26, reduces in the high-energy limit to 
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where 7 = 1/y 1 — /3 2 , and (3 = v/c is the speed of the 
charges, Za and Zb- (7 in the collider frame is related to 
7t in the target frame by jt = 27 2 — 1.) Equation |l]) is 
written in the collider frame, with natural units (c — 1, 
m e — 1, and U = 1), a is the fine-structure constant, and 
a and 7 M are Dirac matrices in the Dirac representation. 
The (S-function form of the interaction is a high-energy 
limit of the exact interaction in a short-range represen- 
tation for the electron's Dirac-spinor, chosen to remove 
the interaction at asymptotic times ]7] ^] . 

For Eq. ([!]) to apply, one assumes first that the ions 
are sufficiently energetic and massive so that the devia- 
tion from straight-line trajectories can be neglected [H. 
Second, one assumes peripheral collisions without nuclear 
interaction. Purely electromagnetic events can be distin- 
guished experimentally from nuclear events by observing 
the full-energy projectile ion after the collision in coinci- 
dence with the produced electrons or positrons One 
also assumes the ions are moving at the speed of light. In- 
deed, for the recent experiments at CERN H, j3 f» 0.99, 
and for future experiments possible at RHIC, j3 « 0.9999. 
Finally, to obtain the (^-function form of the interaction, 
one assumes that 7 ^S> \ f± ± b\, and 26. 

Equation (Q) displays a unique electromagnetic in- 
teraction. The ion with charge Za is moving to the 
right at the speed of light. Its electromagnetic poten- 
tial in the representation chosen here is Lorentz con- 
tracted to the plane transverse to its trajectory, (the 
light front z — t), hence the S(z — t) functional de- 
pendence. Likewise, the ion with charge Zb carries 
with it at the other light front, [z = — t), a plane of 
singular interaction moving to the left at the speed of 
light. Anywhere but on the light fronts, Eq. ([!]) re- 
duces to the free Dirac equation. The Dirac plane waves 
{|x P (f, t)) — e-xjp(-iEpt) exp(ir-p)\up)} which satisfy the 
free Dirac equation are each characterized by the quan- 
tum numbers p = {p, X p , s p }; the momentum p, the sign 
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of the energy E p = (— l) Ap \Jp 2 + 1, and the spin s p = ±. 
Explicit forms for the four four-spinors |u p ) are given in 
Rcf. ||, andp± =p z ±E p . 

The scattering problem for the electron is defined by 
Eq. (|l|) and by plane-wave asymptotic initial and final 
states. One solves for the transition amplitude A% , 

',C0 
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where ^ fc stands for integration over k and summation 
over Afc and Sk- We have obtained an exact, closed- form 
integral representation for the scattering amplitude 
in the following manner 0. First we have observed that, 
as the ions are approaching from infinity, no change oc- 
curs in the region of space between them (\z\ < \t\) until 
t = when the two singular interaction planes collide. 
However, as each ion sweeps through space, it interacts 
with the single plane wave it encounters, resulting in a 
superposition of plane waves after it passes. Each 5- 
function interaction induces a phase-shift discontinuity 
in the wavefunction across each light front J7],||,[l|] . A 
phase shift induced on a plane wave by the passage of a 
single ion is not sufficient to produce pairs, but as the 
two phase-shift planes collide at t = 0, they interfere, 
and pairs are produced as a result. After the collision 
(t > 0), as the ions move apart, the solution in the space 
between them (\z\ < t) is a new superposition of plane 
waves which is determined by the nontrivial boundary 
condition at the light fronts. We have calculated the 

(i) 

transition amplitudes AY by integrating the flux of the 
conserved transition four-current which flows into this 
region across the light fronts. Two terms contribute to 
the amplitude corresponding to the two time orderings 
of the interaction of the electronic wavefunction with the 
two ions. 

The transition amplitudes, A^\ are represented in 
terms of the transverse-momentum transfer distribution 
induced by a single ion, Q b z (K), which contains all the 
dynamics of the ion-electron interaction. When = 
and Xj = 1, A k is an amplitude for a transition from 
the negative continuum to the positive continuum, i.e. an 
amplitude for pair production. We have found Q, 
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The spinor part is o~L(j)±) = (27r) (m^K^ — a z )(a ■ p± + 



\uj) and the momentum-transfer distribution 



is given by the Fourier transform of the phase shift at the 
light front 0, 
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where Z is the charge of the corresponding ion and k 
is the transverse-momentum transfer. The integral over 
£ = n\r±_ — b\ in Eq. should be regularized so as to 
avoid unphysical contributions from large, transverse dis- 
tances, i.e. from £ > K-y. We studied several regulariza- 
tion schemes (for re ^ 0) which all gave M] 
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The exact amplitudes in the infinite 7 limit are ob- 
tained by substituting the result of Eq. (||) in Eq. (||) , 



lim A[ 3) = 



b yi2 a (z A +z B ) T (_ laZA) T (-taZ B ) 



T(+iaZA) T(+iaZ B ) 



x— a 2 Z A Z B / dp ± (p ± - k±) 2 (jp± -J±Y 



a ib-(j±+k±-2p_ l _) 



\px-kxf aZA \p±-Jx\ i2aZB 



Pj_ + 1 - 3+k- 



,—ib-(j± +k± -2p± ) 



Iri 



r aZB \p±~i±\ t2aZA ]} 



(G) 



The branch-point singularities for the intermediate mo- 
mentum p± = k±_ or jj^ are an artifact of using Eq. (|5|) for 
re = 0. (The integral over jfj_ in Eq. (^) has no singulari- 
ties.) We continue the analysis assuming an appropriate 
rcgularization at these points. 

Equation (Q) is nonperturbative, and already includes 
the interaction to all orders in aZ for any value aZ may 
have. Yet, its form is very similar to the high-energy 
limit of results obtained from the two-photon exchange 
diagrams of second-order perturbation theory || , (which 
in the following we simply call the perturbative result). 
For high-energy collisions, there is therefore no reason to 
calculate higher order diagrams. As 7 — > 00, the only 
corrections to second-order perturbation theory calcula- 
tions for free-pair production, including both higher or- 
ders and nonperturbative effects, are the phases in the 
square brackets. For small values of aZ : these phases 
tend to 1 and the perturbative limit is reproduced [0. 
What are the observable nonperturbative effects for fi- 
nite charges? The phase in the square brackets outside 
the integral over p± has no physical implications, but 
the phases in the integrands may substantially alter the 
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physical predictions. We find, for example, that the high- 

/ - \ 2 

energy limit for \Ar'\ differs, in general, from the per- 
turbative result. On the other hand, using Eq. (||) to 

-. , -N 2 

calculate the integrated observable J d(2b)\A^'\ , we get 
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which is identical to the perturbative result. The inte- 
gration over the impact parameter results here in cancel- 
lation of the nonperturbative phases, as J d(2b) exp(i26 • 
p 1 ) = (2tt) 2 5(p > ). We conclude that while some observ- 
ables are sensitive to the nonperturbative phases, other 
observables are not, e.g. because these phases are aver- 
aged to one by an integration. In these cases, observed 
results would agree with the second-order perturbation 
theory calculations, regardless of the size of aZ. 

Would our results apply in an actual experiment, 
where 7 is finite? Equation ([l]) is incorrect for large r± 
(or large b) where it describes an interaction which con- 
tinually increases in strength. Implicit in using Eq. (Q) 
for large, finite 7 is a nontrivial assumption that large, 
transverse distances do not contribute to pair produc- 
tion. In the recent experiments at CERN 7 « 10, 
while in possible future experiments at RHIC and LHC, 
7 w 100 and 7 w 3000, respectively. For these values 
of 7, Eq. (^) is respectively limited to pairs produced 
at transverse distances much smaller than 10, 100, and 
3000 Compton wavelengths away from the ions. This re- 
striction is consistent with an experimental observation 
according to which the average length scale for pair pro- 
duction in relativistic heavy-ion collisions is one Comp- 
ton wavelength B. Yet, being concerned here not with 
averages but with complete distributions, we can not ex- 
clude the possibility of some pairs being produced at 
large, transverse distances from the highly charged ions, 
as long as the transverse-momentum transferred is suf- 
ficiently small. For Eq. (||) to be meaningful, the regu- 
lated integral of Eq. (^) must converge to the expression 
of Eq. (ph for £ such that \r± — b\ <C 7. The case of 
small coupling was previously studied The case of 
large aZ can be considered by the method of stationary 
phase. Expansion of Eq. (^) around the stationary point 
r±_ — b = 2aZn/ k 2 confirms Eq. (||) for this case. The 
procedure is consistent if the stationary point is located 
at small distances from the ion, i.e. if and only if 
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It is interesting to find that Eq. (g) is trivially satisfied 
in two very different limits: in the perturbative limit of 
aZ — > and in the high-energy limit of 7 — > 00. 

Thus, the results which we have first obtained for in- 
finite 7, apply for finite 7 as well. The only restriction 
is of Eq. (||), i.e. that the transverse- momentum trans- 
fer is not too small. For pair production, it is a suffi- 
cient condition to assume that either the initial or final 
(i.e. positron or electron) transverse-momenta are much 
larger than 2aZ/^ where Z is the largest free charge in- 
volved in the collision. The argument goes as follows. 
There are three two-dimensional integration variables in 
Eq. (||) . We first integrate over p± to obtain simple com- 
binations of the Bessel functions of the third kind, Kq and 
K\ . We then use the condition that one of the two trans- 
verse momenta, j± or k± , is much larger than 2aZ/j to 
apply a stationary phase calculation to one of the coor- 
dinate integrations. If, on the other hand, one of the 
charges is screened (a target charge, for example) the in- 
tegral with it converges and there is no need to restrict 
the momentum conjugate to it. The last integral, over 
the other coordinate-integration variable, converges due 
to the Bessel functions which drop exponentially for large 
values of their arguments. Having thus proved that con- 
tributions for the 6-fold integral of Eq. (||) from large, 
transversal coordinates can be neglected, we can make 
the substitution of Eq. (||) and obtain Eq. (^). We re- 
mark that the convergence of the p± integration to the 
Bessel functions occurs only for pair-production ampli- 
tudes for which 1 — j±k T > 0, and is directly related to 
the mass gap between the two continua. It should be 
reconsidered for transitions within the same continuum. 

We now consider the application of our results to the 
discussion of recent, pioneering experiments on pair pro- 
duction performed at CERN's SPS 0. These experi- 
ments measured momentum spectra of positrons emitted 
from pair production in peripheral collisions of 33-TeV 
Pb 82 + ions (7T = 168) and 6.4 TcV S 16+ ions ( 7T = 212) 
with various targets (i.e. (CH^xj Al, Pd, and Au). The 
charge dependence of the positron yield was reported 
with excellent precision. The target-charge dependence 
for the sulfur projectile is Z^ 99±0 ' 02 , and for the lead pro- 
jectile is Z^ 03 ^ 03 ; both within - 1% agreement with the 
prediction of perturbation theory. The projectile-charge 
dependence was observed to be Zp 0±01 , also in very good 
agreement with perturbation theory. The positron mo- 
mentum distributions for sulfur and lead projectiles are 
compared by scaling each spectrum by Z|>, and by scal- 
ing the sulfur data from 7t = 212 to 7t = 168 as ln 3 (7-r), 
as predicted by perturbation theory. The scaled distribu- 
tions are observed to be approximately the same, and to 
agree reasonable well with two-photon perturbation the- 
ory (see discussion in |Q), except for enhancements for 
the lead projectile at very low (< 2 MeV/c) and high (be- 
tween 8 MeV/c and 12 MeV/c) momentum. The authors 
of Ref. P| note that the variation of the scaled momen- 
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turn distribution with the projectile charge, and not the 
target charge, is unexplained. 

The observed Z^Z" 2 ? charge dependence of the single- 
positron yields, even for very large charges, is consis- 
tent with the charge dependence we have obtained for 
the nonperturbative, high-energy limit (see Eq. (Q)). It 
agrees with perturbation theory but is not a perturbative 
effect. Nonperturbative phases in the exact amplitudes 
make them different from second-order perturbation the- 
ory results, but these phases cancel for calculations of 
total cross sections. Our theoretical prediction of a Z 2 
dependence of the total cross section in the high-energy 
limit implies that multiple-pair production in very high- 
energy collisions cannot be inferred from a measurement 
of the charge dependence of the total positron yield . 
We suggest that the two regions of excess cross section 
observed in the experiment have a common origin: an 
enhancement over perturbation theory for small values 
of the transverse-momentum transfer, for which Eq. (Q) 
diverges. We found that agreement with the perturba- 
tive result is restricted by Eq. (||) . Assuming that in the 
collider frame pairs are produced isotropically, this re- 
striction, formulated for the transverse momentum, may 
translate to a restriction on the total positron momen- 
tum: j± ~ j z ~ j >> 2Zpa/"f. Taking >> to be a 
factor of 10, we then predict, in very good agreement 
with the observed scaled spectra, that the perturbative 
result for the positron momentum distributions is valid 
for j > 0.4 MeV/c for the sulfur data, and j > 2 MeV/c 
for the lead data. The excess cross section observed at 
high momentum for the lead projectile is consistent with 
an enhancement in the cross section at low transverse 
momentum in the collider frame after a relativistic trans- 
formation to the target frame is applied. The absence of 
an observed target-charge dependence for the scaled dis- 
tributions, is most likely attributable to screening by the 
atomic electrons ifTlj] . 

In conclusion, we have shown that the exact, non- 
perturbative solution and the two-photon exchange di- 
agrams of second-order perturbation theory give exactly 
the same results for free-pair production yields integrated 
over the impact parameter, as long as the transverse- 
momenta transferred from the ions to the electron are 
larger than 2aZ/j. The leading-order perturbative cal- 
culations for this observable are therefore exact not only 
at the perturbative limit of aZ -C 1 but also in the high- 
energy limit of 7 3> 1. This explains recent experimen- 
tal results according to which production rates scale as 
ZpZ|,, even for large charges. New nonperturbative ef- 
fects could be detected by measuring observables differ- 
ent from the integrated, inclusive production rate that 
was measured in these experiments. The exact ampli- 
tudes of Eq. (JsJ) include nonperturbative phases which 
may have an observable effect, e.g. if one does not in- 
tegrate over the impact parameter 26. We expect these 
phases to strongly influence the theoretical predictions 



for correlations and multiple-pair production. Several is- 
sues deserve further study. These include pair production 
at large, transverse distances from the ions and bound- 
free production for which one should obtain the solution 
on the light fronts themselves. For other calculations, 
Eq. (^) as well as the physical picture that led to it, are 
likely to become useful theoretical tools |l2| . 
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